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We derive a quantum master equation from first principles to describe friction in one dimensional, 
collisional Brownian motion. We are the first to avoid an ill-defined square of the Dirac delta function 
by using localized wave packets rather than plane waves. Solving the Schrodinger equation for two 
colliding particles, we discover that the previously found position diffusion is not a physical process, 
but an artifact of the approximation of a coarse grained time scale. 
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More than a hundred years ago, Einstein and Smolu- 
chowski [T] derived a successful classical description of 
Brownian motion, the erratic movement of a particle af- 
fected by ambient gas molecules. In the last few decades, 
the field of quantum Brownian motion (QBM) has be- 
come established [2J, [3]. Not only is such a quantum 
description of friction and diffusion desirable, but QBM 
also sheds light on the non-classical process of decoher- 
ence which is believed to be of importance in the quantum 
classical transition. 

The field of QBM is still the source of open questions. 
Of particular concern is deriving a valid master equation 
(ME) for the reduced density operator of the Brownian 
particle in the Markov limit, in which case a ME of Lind- 
blad form is to be expected. What should be noted here 
is that all such MEs include position diffusion. This pro- 
cess is impossible in classical dynamics (although the po- 
sition does diffuse in the long term due to randomness of 
the momentum) as it derives from position jumps which 
are at odds with continuous trajectories of classical par- 
ticles. Therefore it is currently believed that position 
diffusion arises as a quantum effect accompanying colli- 
sional quantum friction 0J. In this letter we show the 
contrary: We find that position diffusion is not a real 
physical process in collisional QBM, but results from us- 
ing a coarse grained time scale on which collisions ap- 
pear instantaneous. Unproblematic in classical dynam- 
ics where collision times of hard-core particles are indeed 
short, the assumption of instantaneous collisions has to 
be used with care in quantum dynamics, where collision 
times depend on the widths of the colliding wave packets. 

There are several "non-collisional" models of QBM in 
the literature (see e.g. [S]-[7]). Here we concentrate 
on collisional models which are the closest in spirit to 
the original Einstein approach, and which can be tested 
in experimental setups [5]. Important contributions in 
the literature include [9]- [12] and can be reviewed in [4]. 
Mostly, the analysis is based on using scattering theory 
to describe the effects of a single collision. 

To point out some problems which usually occur in 
collisional QBM when using scattering theory, we briefly 
review Hornberger's contribution |12j as it seems to be 



the most complete one. Due to momentum conservation 
a two particle collision reduces to a one particle problem, 
where the transition operator T with momentum matrix 
elements related to the scattering amplitude (p'|T|p) = 
5(p a — p 2 )/(p',p)/(7rfi-) is of importance. In particular 
it appears in an expression of the form 
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(p' + q|T|p + q)(p'-q|Tt|p-q) 
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where f2 — > oo is a box-normalization volume. For q = 
this term contains an ill-defined square of the Dirac delta 
function and a physically motivated replacement rule 
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is used. However, for q ^ Eq. (JT|) is well defined. In 
fact, using standard relations for Dirac delta functions, 



X = 



|UtP-^((P- P 'W 

x f(p' + q, p + q)/* (p' - q, p - q) 



(3) 



is found, where q = |q| and p q is p projected on q. Nev- 
ertheless, this result Eq. (|3| is not without its problems as 
the off-diagonal matrix elements (p'|p|p) of the Brown- 
ian particle's density matrix would be found to disappear 
in a single collision due to f2 — * oo, no matter how small 
p — p'. This seemingly strange observation comes with 
the fact that the collision time with a gas particle in a 
non-localized momentum eigenstate is infinite. 

To construct a Markovian (memory less) ME, Horn- 
berger assumed instantaneous collisions. But within this 
assumption, the complete vanishing of coherences is not 
acceptable any more. To get around this problem, Horn- 
berger substituted the square root of the replacement rule 
Eq. @ for (p'±q|T|p±q) in Eq. Q to effectively bring 
back the otherwise lost coherences of momentum states. 

The replacement rule Q has to be used with caution 
in the first place, but its application to a well defined 
expression is even more questionable. Whether decoher- 
ence in momentum bases and the closely related position 
diffusion are correctly described by approaches of this 
kind is surely not certain. 
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To resolve this matter we study a single collision in 
terms of localized gas states. Then a collision time 
can be precisely defined, and the low-density and high- 
temperature (LDHT) limit will be quantified by 'colli- 
sion time* x ' collision rate' <C 1. To examine a collision 
event time resolved, and to investigate under which con- 
ditions a complete collision of two wave packets occurs, 
we analytically solve the two-particle Schrodinger equa- 
tion. This is in contrast to scattering calculations which 
postulate a complete collision. 

To simplify the problem, we consider one-dimensional 
QBM and assume a hard core interaction potential be- 
tween gas and Brownian particles. Our notations are as 
follows. Position and momentum operators are denoted 
with a hat, whereas wave function variables are primed. 
The index g refers to gas particles. 

The two particle Hamiltonian 



H 



A. 

2m„ 



+ h ad(x - X g ), 

2m 



(4) 



forbids any tunneling of the gas particle through the 
Brownian one. The two orthogonal sets of energy eigen- 
states in position space are 

i>t h (x' g ,x') = e i ^' +ax 'shm(k(x' g -x% (5) 
= e^'+^Osin^l^-x'l), (6) 

with k € R, k € K+, a — m g /m, and energy E(k, k) — 
(1 + a)(k 2 + ak 2 )h 2 /{2m g ). 

The collisional analysis is based on noting that the 
density operator of a thermal ideal gas particle at tem- 
perature T can be written as convex decomposition of 
Gaussian states 



Thus we have to construct an initial two-particle state, 
expanded in terms of energy eigenstates Eq. ([5| and (|6| , 
which satisfies the boundary condition ip(t = 0,x g ,x = 
x g ) = required by the hard core interaction, and which 
approaches Eq. (|9| in the limit Eq. (10). As the proce- 
dure is quite lengthy, we only present the result and refer 
to our more detailed paper [2]. It may be checked that 
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has indeed the required properties. Without loss of gen- 
erality we used the reference frame of centre of mass 
(i.e. p = —p g and x = —ax g ) as well as x g < x and 
p g > p. Furthermore, we related the widths aa 2 = a 2 
of the Gaussian wave packets to the relative mass of the 
colliding particles. This turns out to be the crucial re- 
quirement for the two particles not being entangled after 
the collision (see also [15]), and \x,p) a can be seen as 
the pointer bases of a measurement performed by a gas 
particle in the state \x g ,p g ) [H] , 

The time evolution is no w e asily achieved by multi- 
plying the integrand of Eq. |ll| with e ~^E{k,k)t/h _ After 
straight forward integration we find 
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dpglI (Jg (Pg)\Xg,Pg} iT (Xg,Pg\, (7) tf)(t,x'x') 



where L is a normalization length and where \x,p) a de- 
notes a Gaussian wave packet with position variance a g : 



{x'\x,p) an = 



- ) —ixp/2h 
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e ix'p/h e -(x-x') 2 /2a 2 g ^ 



Hornberger and Sipe [T3] showed that Eq. ([7| is valid 
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We therefore undertake this analysis by considering an 
initial state that approaches the product state 

\Xg,p g ) O \x,p) (9) 

in the limit of large separation, i.e. we assume vanishing 
overlap of the two initial wave packets 



\Xg ~ X\ » 



(10) 
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for x' g < x' and ip(t, x' g , x') — otherwise. At time t = 
the first term in the curly brackets is the dominant one. If 



relation (10 1 holds, then we can neglect the second term 



and Eq. (12 1 becomes the position representation of the 
desired initial state (J9j> . If the particles' relative velocity 
is large compared to their velocity uncertainty 



\Pg\ » 



(13) 
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then we can neglect the first term sufficiently long after 
the collision. In this case Eq. (12) is the position repre- 



sentation of the remarkable simple product state 

Ug(t) \-Xg,- Pg )®U(t) \~X,-p) , 



(14) 



where U(t) is the free particle evolution operator. Using 
Eq. (Il2|) we can now define the collision time 



tr. = 



1 



(15) 



as the time it takes from the first term in the curly brack- 
ets to dominate, to the second term to dominate. 

To examine the behavior of the Brownian particle 
during a collision in detail we trace out the gas parti- 
cle from Eq. ( fl2| ). Its position probability distribution 

P (t,X') = jdX^ ) can be obtained analyti- 

cally in terms of the error function and does not show 
any position jumps. An example is plotted in Fig. [TJa), 
and reminds of a classical collision. As Gaussian states 
build an over complete set of states, this result shows 
that no position jumps, and hence no position diffusion 
can occur in collisional QBM. Using the Fourier trans- 
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FIG. 1: (color online) Position (a) and momentum (b) prob- 
ability distributions for the Brownian particle during a col- 
lision. While the position probability distribution "flows", 
the momentum distribution "jumps" from the initial value to 
the final one. This disproves position diffusion in collisional 
QBM. Parameters are: x = 10, p = —2, m = 1, a = 0.3, 
h = 1, and cr = 4. 



formation of Eq. (12 1 with respect to x', we find the 



momentum distribution of the Brownian particle which 
shows the expected momentum jump (see Fig. [ljb)). 

In a general reference frame the collision event with 
Gaussian states can be written as 

\Xg,Pg)®\x,p) " CattCrmS > Ug(t)\Xg,Pg)®U(t)\X,p) (16) 

where positions and momenta after the collision relate 
with the initial values in the classical way |17j : 



2x — (1 — a)x 
l + c 



- = 2ax 9 + (1 ~ a ) x 



_ 2ap - (1 - a) Pg ^_ 2 Pg + (l-a)p 



The transformation of a general Brownian particle's den- 
sity matrix due to a collision with a gas particle in the 
state \xg,p g ) can be written in terms of Kraus operators 

p(t) = JJdMpU(t)K Xg , Pg (x,p) P (0)Kl. Pg (x,p)UHt). 

(19) 

In |14j we will present a derivation of the Kraus operators 
K XgPg {x,p) = d(^-{x 9 -x), j^—(p g -ap^j ^tc(x,p) 



where D(x,p) 
operator and 

n(x,p) 
w(x,p) 



_ p i(px-xp)/H 



is the Glauber displacement 



dxdp 
2nh 
2a 



w(x,p) \x + x,p + p) (x + x,p + p 
2a 



nh(l — a)'- 



exp 



h 2 



are effect operators corresponding to a measurement 
which is indirectly performed by the gas particle [14] . 
The transformation (19 1 may be confirmed by applica- 



tion to Gaussian states and comparison to Eq. (16 1 



Next we determine whether during a time interval 5 
the Brownian particle actually collides with a gas particle 
which is in the state \x gi p g ). For this purpose we assume 



S > t c 



(20) 



so that we can neglect partial or uncompleted collisions. 
Then we find that the two particles collide if (x,p) is 
within the phase space region S(x g ,p g ) = {(x,p)\0 < 
— , x Zg ; — < S\. Therefore we have to project the Brow- 
nian particle's density operator onto S(x g ,p g ) before ap- 
plying the transformation Eq. (19). Of course, phase 



space projections do not exist in quantum mechanics, but 
as long as ( 20 1 is satisfied we can use the approximate 



projection operator 

T 8(Xg,Pg) = 



S(Xg, Pg ) 



dxdp 
2nh 



x,p)(x,p\ 



(21) 



By multiplying this operator with the probability of 
finding a gas particle in the state \x g ,p g ), we define the 
collision probability operator 



P s(Xg,Pg) = ngP(Pg)T S (Xg,P g ), 



(22) 



where n g is the gas number density. In fact, the expecta- 
tion value Tr[pPs(xg,p g )] gives the probability of a colli- 
sion of the Brownian particle and a gas particle described 
by \x g ,Pg). We also define Pg — J J dx g dpgPs(x g ,Pg) cor- 
responding to the total collision probability. 

To find the density operator for the Brownian particle 
at time <5, quantum trajectory theory [16] suggests to 
multiply each Kraus operator by the square root of the 
probability operator. We use the low density assumption 



Tr(pPs) « 1, 



(23) 
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to neglect two or more collisions during 6 and find 



p(8) = U(6) 



dxgdp g dxdp K XgtPg (x,p) 



x^P s (xg,p g )p(0)^Ps(x g ,p g )Kl gPg (x,p) 
+ ^1-P 5 p(0)^l-P s } tf(6). 



(24) 



The usual approach to a ME is to expand all operators 
in 6 and to take the limit 8 — > 0. However, this limit is 
not allowed in Tg(x g ,p g ) as then uncompleted collisions 
were not negligible. A finite 5 will ultimately lead to 
errors which we minimize by performing the transition 
to continuity in the interaction picture and to second 
order in 8 



Pint(8) - Pmt{0) _ Pint ( f + dt) - p mt (f ) 



6 dt 
Transforming back to the Schrodinger picture we find 



(25) 



dp(t) 
dt 



~[H,p(t)]-^[Rp(t) + p(t)R] 



dx g dp g dxdp U(\)K x ^ Pg (x,p) 



jRs{x g ,p g )rt {l)p(t)U(l ) ^Rs{x giPg ) 



(26) 



x Kt g>Pg (x,p)Um) 



where R = Pg/8 and R${x g ,p g ) = Ps(x g ,p g )/8 are rate 
operators corresponding to collision rates. This master 
equation of Lindblad form is our main result. It describes 
collisional friction as well as dccohercncc for general mass 
ratios and velocities, and allows a clear interpretation in 
terms of quantum trajectories. 



An interesting feature of Eq. ( 26 1 is the dependence on 
d and <j g . These parameters persist in the ME because 
they are bound from above and below by the relations 



(10 1 and (13), and (20 1 and (23 1, respectively and there- 
fore the limits 8 — *■ and a g — > oo are not allowed. 
However, we show in [T3] that the effects of 8 and a g on 
the dynamics vanish in the LDHT limit 



V2(l + a)^^= « 1, 



(27) 



which is also the crucial condition to simultaneously sat- 
isfy relations (10 1, ( [T3| , (20 1, and (23 1 for most thermal 
gas momenta. Our ME is valid on a coarse grained time 
scale 8 > (1 + a)h/(k B T). 

We apply our ME to derive equations of motion for 
expectation values of position and momentum. In the 
limit of a slow Brownian particle \p/m\ <C ^Jk B T /m g 
(the generalization to fast Brownian particles is found in 
14J and includes non- linear friction) we find 



d{x) 
dt 



(P>. 
m 



dt 

d(x 2 ) 

dt 

d({x,p}) 
dt 

d(P 2 ) 
dt 



-m, 

({x,p}) n 



gas 



3^ 



2k B T 



3/2 



Hp 2 ) 



-imp})- 



= 2f[mk B T~(p 2 )] 



(29) 
(30) 

(31) 
(32) 



where / = 4n g y /, 2m g k B T / \y/itm) is the friction con- 
stant. The second term in Eq. p0] ) represents position 
diffusion. However, as shown in our single collision cal- 
culation this is not a true physical process. It rather 
results from neglecting partial or uncompleted collisions, 
which is necessary for any Markovian description of col- 
lisional QBM. This is confirmed by the observation that 
this term is small in the range of validity of our ME. 

Using only first principles we derived a quantum ME 
of Lindblad form for collisional QBM in one dimension 
in the LDHT limit. We showed that the counter intu- 
itive position diffusion found in earlier work results from 
approximations necessary to derive Markovian dynam- 
ics, but is not a real physical process. The application 
to expectation values, where we recovered the expected 
classical equations served as a test of our ME. In a fu- 
ture paper we will use our ME to study the important 
issue of collisional decoherence which is now in reach of 
quantitative experiments [5]. 
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